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First Lattice Study of the iV-Pn( 1440) Transition Form Factors 
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Experiments at Jefferson Laboratory, MIT-Bates, LEGS, Mainz, Bonn, GRAAL, and Spring-8 
offer new opportunities to understand in detail how nucleon resonance (N*) properties emerge from 
the nonperturbative aspects of QCD. Preliminary data from CLAS collaboration, which cover a 
large range of photon virtuality Q 2 show interesting behavior with respect to Q 2 dependence: in the 
region Q 2 < 1.5 GeV 2 , both the transverse amplitude, Ai/2(Q 2 ), and the longitudinal amplitude, 
Si/2(Q 2 ), decrease rapidly. In this work, we attempt to use first-principles lattice QCD (for the first 
time) to provide a model-independent study of the Roper-nucleon transition form factor. 

PACS numbers: 11.15.Ha, 12.38.Gc 13.40.Gp 13.40.Hq 14.20.Gk 



Introduction: Lattice QCD has successfully provided 
many experimental quantities from first-principles calcu- 
lation; however, its success has mostly been restricted to 
measurements of ground-state quantities. Lattice mea- 
surements of excited states could contribute, for exam- 
ple, to hadron spectroscopy, where there are many poorly 
known states which require theoretical input to be identi- 
fied. At the EBAC at Jefferson Lab, dynamical reaction 
models have been developed to interpret extracted N* 
parameters in terms of QCD0,[3. 

Among these excited nucleon states, the nature of the 
Roper resonance, AT(1440) or N' , has been the subject 
of interest since its discovery in the 1960s. It is quite 
surprising that the nucleon's excited-state mass is lower 
than its opposite-parity partner, a phenomenon never ob- 
served in meson systems. There are several interpreta- 
tions of the Roper state, for example, as the hybrid state 
that couples predominantly to QCD currents with some 
gluonic contribution Q or as a five-quark (meson-baryon) 
stateQ. Some earlier quenched lattice QCD calculations, 
e.g. Refs. [1, @, 0, S, , find a spectrum inverted with 
respect to experiment, with N' heavier than the opposite- 
parity state Sn. However, in the study of Ref. [ll|, where 
larger lattice box and lighter pion masses are used, their 
findings of indicate a rapid crossover of the first positive- 
and negative-parity excited nucleon states close to the 
chiral limit. The lattice study has not ruled out that the 
roper is the first-excited state of nucleon and this is our 
assumption in this work. 

The new data provided from experiments at Jefferson 
Laboratory, MIT-Bates, LEGS, Mainz, Bonn, GRAAL, 
and Spring-8 offer a new opportunity to understand in 
detail how nucleon resonance (N*) properties emerge 
from the nonperturbative aspects of QCD. For exam- 
ple, the extraction of the — ► N* transition form 
factors could help us to understand the dynamical ori- 
gins of the confinement of constituent quarks and the 
associated meson cloud. There are various QCD-based 
hadron models, such as the well developed constituent 
quark model 12, [l3| and the covariant model based on 
Dyson-Schwinger Equations 



Understanding of the true nature of the Roper may be 
most easily gained by studying its structure and form 
factors, such as the nucleon-Roper transition. CLAS 
collaboration 15, H, [l?], EH has studied such transitions 
induced by electron scattering over a large range of in- 
termediate photon virtuality Q 2 . In the region Q 2 < 
1.5 GeV 2 , both the transverse amplitude, Ai/ 2 (Q 2 ), and 
the longitudinal amplitude, S\/2(Q ), drop rapidly in 
magnitude. This is well described in relativistic quark 
models with light-cone dynamics, and the sign is consis- 
tent with the non-relativistic version. However, in the 
low-Q 2 region, A 1 / 2 (Q 2 ) becomes negative; this is not 
understood within constituent quark models and requires 
inclusion of meson degrees of freedom. At the EBAC cen- 
ter, dynamical reaction models have been developed to 
interpret extracted N* parameters in terms of QCD 0,0]. 
A model-independent study of these quantities from lat- 
tice QCD will serve as valuable help to phenomenologists 
in analyzing experimental data and will provide better 
theoretical ground for understanding low-Q 2 physics. 

The AT(1440) — * N Form Factors: From Lorentz sym- 
metry, we expect the matrix element composed of any 
nucleon states, Ni and N 2 , to have the following general 
form: 
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where the equation of motion is used to simplify —5^7^ = 
M N2 -M Nl . 

To compare with the experimental results, we connect 
the experimentally measured helicity amplitudes A 1 / 2 
and S1/2 to the transition form factors F* 2 through 

A 1/2 (Q 2 ) = k A (Q 2 )G M (Q 2 ) 

S 1/2 (Q 2 ) = k s (Q 2 )G E (Q 2 ) (2) 

with 
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k s (Q 2 ) ^ MQ 2 ) f % + ^ a ^ Q2 + ~ Mn ^ 

x yjQ 2 + (M Nl +M N2 ) 2 , (4) 
the magnetic and electric transition form factors 



(0\(x N y\p, s) — Zu s (p) with (with nucleon interpolat- 
ing field X )i the spinors satisfy 
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g m (Q 2 ) = f*(q 2 ) + f;(q 2 ) 

G e (Q 2 ) = F*(Q 2 ) 
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Traditionally, one is only interested in the ground state; 
thus the smearing parameters are chosen to overlap as 
little as possible with excited states. We will use var- 
ious smearings to create correlators having measurable 
overlap with the Roper excited state. 
Similarly, the three-point function is 
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and QED running coupling a ~ 1/137. Using these defi- 
nitions, we can reconstruct F* 2 from experimental values 
of helicity amplitudes. 

Lattice Setup: In this exploratory study, we will 
use an anisotropic lattice; that is, a lattice where the 
temporal lattice spacing is finer than the spatial ones. 
It has been demonstrated in the past that for certain 
calculations, such as glueballs[l9| and multiple excited- 
state masses [20|, there are great advantages to using 
anisotropic over isotropic lattices due to the finer lattice 
time spacing, even when the fundamental constituents 
are not heavy. 

We perform our calculations on quenched 1 16 3 x 64 
lattices with anisotropy £ = 3 (i.e. a s = 3a t ), using Wil- 
son gauge action with f3 = 6.1 and stout-link smeared 21 
Sheikholeslami-Wohlert (SW) fermions[22j (with smear- 
ing parameters {p,n p } — {0.22,2}). The parameter v is Zy is the vector-current renormalization constant and FF 



xe -(t/-t)K(p/) e -(*-*0^(F l ) x ME's (8) 
where n and n' index energy states and 

ME's = Y^ AE n'{Pf)E n {Pi)T a pu n >(p f ,s') 

s,s' 

x (N n >(p f ,s') \V^\N n {p t) s))u n {p tl s) a 
= Tr [T ■ (E n ,(p f rf - *7 ■ py + M n /) 

•FF's- ( J B„(p l ) 7 *-*7'P l + Af„)] , (9) 

where the projection used is T = |(1 + 74)(1 + 27573), 



nonperturbatively tuned using the meson dispersion rela- 
tion, and the clover coefficients are set to their tadpole- 
improved values. The inverse spatial lattice spacing is 
about 2 GeV, as determined by the static-quark poten- 
tial, and the simulated pion mass is about 720 MeV. In 
total, we use 200 configurations of anisotropic lattices. 

In this work, we will use the variational method and 
simultaneous fitting on two- and three-point Green func- 
tions to extract N-Pu transition form factors. The nu- 
cleon two-point correlators measured on the lattice are 



E n + M„ 
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where A and B index over the different smearing pa- 
rameters and n indexes over the basis of nucleon energy 
eigenstates; the states are defined to be normalized as 



denotes the form factors. The vector current in Eq. [5] is 
0(a) on-shell improved with the improved coefficient set 
to tree- level value 2 . Zy is calculated nonperturbatively 
from the isovector vector charge Zy = l/g v ~ d - Note 
that when one calculates the three-point Green function 
in full QCD, there are two possible contraction topolo- 
gies: "connected" and "disconnected" diagrams, when 
the vector current vertex appears on a vacuum bubble. 
In this work, only "connected" quantities are included. 

Numerical Results: We apply a smearing function 
with gauge-invariant Gaussian form to improve the over- 
lap of fermion operators with states of interest. To ob- 
tain matrix elements involving more than ground states, 
we need to be careful not to over-smear the fermions; 
a widely smeared fermion will greatly suppress excited- 
state signals. (This goes against the standard prac- 
tice used in the nucleon matrix elements where one 
increases the smearing parameters to suppress excited 
states.) We use three Gaussian smearing parameters: 



A "quenched approximation" means the effective sea-quark mass 
is infinite. The quenched QCD (QQCD) is confining, asymp- 
totically free, shows spontaneous chiral symmetry breaking and 
differs from full QCD only in the relative weighting of the gauge 
configurations, it is reasonable to use quenched simulations to 
test new lattice techniques. QQCD is very useful for under- 
standing and controlling sources of errors (except for that due to 
the quenched approximation) for a new methodology before it is 
extended to a much more expensive full-QCD calculation. 



2 With smeared fermion actions, it has been seen on three-flavor 
anisotropic lattices with tree-level tadpole-improved coefficients 
in the fermion action that the nonperturbative coefficient con- 
ditions are automatically satisfied 23] . Similar behavior has also 
been observed in a quenched studv|24ll. where the nonpertur- 
bative coefficients or renormalization constants in a smeared 
fermion action differed from tree-level values by a few percent. 
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(7 e {0.5, 2.5, 4.5}; the largest of these has excellent over- 
lap with just the ground state, while the other two in- 
clude substantial higher-state contributions. For both 
two-point and three-point correlators, we calculate all 9 
possible source-sink smearing combinations. 

From the two-point functions, we wish to extract the 
energies of at least the nucleon and Roper (P's) and the 
matrix elements of our lattice operators between the vac- 
uum and the baryon states (2Ts). We apply the varia- 
tional method[25} to extract the best principal correlators 
corresponding to pure energy eigenstates from our matrix 
of correlators. A simple exponential fit to these yields 
M N = 1.48(2) GeV and M R = 2.53(8) GeV. The Z's 
are then determined from the eigenvectors of the princi- 
pal correlators; we select the eigenvectors from the time 
tz that minimizes the discrepancy between the correla- 
tors reconstructed from our Z's and P's and the origi- 
nal three-point data. We find the mass of the negative- 
parity partner Sn using the same method: Ms 1 i = 
2.40(10) GeV. This mass is lower than the Roper mass 
due to the high quark masses used, a result consistent 
with previous lattice calculations; it is expected that the 
masses will cross when lower quark masses are used. 

We calculate three-point functions with source and 
sink locations at 15 and 48 respectively. (This gives a 
source-sink separation about 1.1 fm.) The final state al- 
ways has zero momentum for convenience, while the mo- 
mentum of the initial state varies. We can get the nucleon 
form factors, F^ N and F^ N , from the large Gaussian 
smearing (er = 4.5) runs using a ratio approach[26j. We 
use the ground-state term of the fitting form in Eq. [8] on 
smearing parameter a = 4.5 three-point correlators. The 
fit range is adjusted so that the fitted results are consis- 
tent under small perturbations to the range. We obtain 
nucleon-nucleon form factors F\ and F% consistent with 
those derived from the ratio method. Note that since 
our pion mass is far heavier than the physical one, we do 
not expect to see very good agreement with experiment. 
As the pion mass approaches lighter values, the lattice 
data will trend toward the experimental values; for more 
details, see the recent review in Ref. 
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The ratio method will not extract matrix elements be- 
yond the ground state, but using the Z's and E's derived 
from our analysis of the two-point functions, we can ex- 
tract excited matrix elements from the three-point func- 
tions by fitting to the form given in Eq. [8j We want to 
keep terms in Eq. [8J having n and n' running from the 
ground to the first-excited state; thus, there would be 
four matrix elements in the minimum expansion, which 
will be free parameters in our fit. We increase the num- 
ber of three-point correlators, first using just the diag- 
onal correlators where the smearing is the same at the 
source and sink, and then using all 9 smearing combina- 
tions. The nucleon-nucleon matrix elements are verified 
against the ratio method, and we check that the tran- 
sition matrix elements are consistent between different 
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FIG. 1: Proton-Roper form factors FTn obtained from CLAS 
experiments 

[IS [M[l2l, El and PDGQ number (circles) and 
our fitting method (squares, diamonds) 



sets. In the rest of this work, we show the results from 
the full 9-correlator simultaneous fits. 

On the lattice, we can obtain both the form factors 
related to the Roper decay Pn — > -fN, and the one that 
is related to photoproduction, j*N — >• Pn. A Roper at 
rest has a P-wave decay into a pion and a nucleon and 
an S- wave decay into two pions and a nucleon via Sii7r. 
Since we do not wish to consider the complicated case 
that occurs when two-particle states may be present, we 
must avoid kinematical situations in which decays can oc- 
cur. In our simulation, these are suppressed by the high 
quark mass and discretization of momentum. However, 
the lowest-momentum Roper {Er ~ 2.7 GeV) can decay 
into a lowest-momentum pion (E v w 1.0 GeV) and a nu- 
cleon at rest {Mm m 1.5 GeV). Since we cannot untangle 
the Roper from the two-particle state here, we drop this 
data point. All other decays are forbidden. 

In Figure [I] we show the transition form factors F* ^ 



derived from experimental helicity anrplitudes[15l. Il6l. 1171 
[l8| and those from our numbers. In Figure [2l we show 
the neutron-Roper transition form factors with an exper- 
imental point. Note that since our nucleon and Roper 
masses are much higher than the physical ones, we are 
in the time-like region when we use the matrix element 
(Pn|V^|iV) to construct the transition form factors. As 
we decrease the pion mass, we will enter the space-like 
region and this matrix element will be helpful in giving 
us different Q 2 points. Our calculation seems to be quite 
different from the experimental values. This is somewhat 
expected, since our pion mass is much heavier than the 
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FIG. 2: Neutron-Roper form factors F* 2 obtained from 
PDG[28|] number (circles) and our fitting method (squares, 
diamonds) 



physical value. Even in unquenched calculations, the nu- 
cleon form factors do not agree with experimental values 
with pion mass as low as 300 MeV[27|. We may expect to 
see the lattice data approaching experiment as the pion 
mass used is decreased in the future calculations. 

Conclusion and Outlook: This exploratory study 
demonstrates that the nucleon- Roper transition form fac- 
tors can be measured from a first-principles lattice cal- 
culation. Using the fitting approach with appropriately 
chosen operator smearing, we not only improve the signal 
in the nucleon-nucleon form factors (especially at large 
momenta), but also successfully extract the nucleon- 
Roper. We may vary the projector used in the three- 
point function to further improve the signal. In the fu- 
ture, since the pion mass in our simulation is very heavy 
at 720 MeV, we would like to consider lighter pion masses 
as well as starting work on unquenched anisotropic lat- 
tices. 
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